Microscopic description of oblate-prolate shape mixing in proton-rich Se isotopes 
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The oblate-prolate shape coexisting/mixing phenomena in proton-rich 68'''0,72gg ^^.^ jjiyestigated 
by means of the adiabatic self-consistent collective coordinate (ASCC) method. The one-dimensional 
collective path and the collective Hamiltonian describing the large-amplitude shape vibration are 
derived in a fully microscopic way. The excitation spectra, i3(E2) and spectroscopic quadrupole 
moments are calculated by requantizing the collective Hamiltonian and solving the collective 
Schrodinger equation. The basic properties of the coexisting two rotational bands in low-lying 
states of these nuclei are well reproduced. It is found that the oblate-prolate shape mixing becomes 
weak as the rotational angular momentum increases. We point out that the rotational energy plays 
a crucial role in causing the localization of the collective wave function in the 7) deformation 
space. 

PACS numbers: 21.60.-n, 21.60.Ev, 21.10.Re, 27.50.-|-e 



I. INTRODUCTION 

Atomic nuclei exhibit various intrinsic shapes in their 
ground and excited states. Coexistence of different 
shapes in one nucleus is widely observed all over the 
nuclear chart Among varieties of such phenomena, 
much attention has been paid on proton-rich N = Z 
nuclei in the A 7Q region where very rich shape coexis- 
tence phenomena are seen. In this region, dramatic com- 
petition of different shapes occurs due to shell-structure 
effects; the oblately deformed shell gaps at or Z = 34 
and 36, the prolately deformed shell gaps at 34 and 38, 
and the spherical shell gap at 40 0. 

The N = Z nucleus ^^Se is a particularly interest- 
ing nucleus, because it has deformed shell gaps both at 
the oblate and prolate shapes. For this nucleus, mean- 
field calculations predict the oblate-prolate shape coex- 
istence [1, 0, Q. In experiment, two rotational bands 
were observed in low-energy excitations, and the ground 
and excited bands were interpreted to have the oblate 
and prolate deformations, respectively jH, 0|- For ™Se 
and "^^Se, B(E2) values for the 2"[,A'i[,Q"[ states in the 
ground bands were obtained by a recent life-time mea- 
surement ^ . These data indicate gradual change of their 
characters from oblate to prolate with increasing angu- 
lar momentum; it occurs in lower angular momentum in 
''^Se compared to ''"Se. The data also suggest consid- 
erable mixing of the oblate and prolate shapes in these 
low-lying states. We also note that candidates for the ex- 
cited O2 states have been known for a l ong time at 2011 
keV in ™Se Q and at 937 keV in ^^Se 

Since shape mixing is caused by large-amplitude col- 
lective motion connecting different shapes, its theoreti- 
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cal description is beyond the static mean-field approx- 
imation or the small-amplitude fluctuation about equi- 
librium shapes. A difficulty in theoretical description of 
the shape coexisting/mixing phenomena is that various 
kinds of microscopic configurations associated with dif- 
ferent shapes participate in them and thus quite a large 
number of particle-hole degrees of freedom are involved 
in the large-amplitude collective dynamics. Therefore, 
microscopic description of shape mixing is a challenging 
subject in nuclear structure theory. 

Theoretical investigations on the shape coexist- 
ing/mixing phenomena may be divided into two cate- 
gories: i.e., time-independent and time-dependent ap- 
proaches. For the former, we can refer, e.g., t he p ro- 
jected shell model [TT|, the lar ge-s cale shell model [l2l.ll3j| 
the interacting boson model [IJI calculations for ^*Se, 
and the excited- vampir variational calculation for ^^'^°Se 
(Tsl [T6| . For neighboring isotopes ^^^^^Kr, a detailed 
study based on the number and angular-momentum 
project ed g enerator-coordinate method was recently re- 
ported [l7| . 

A well-known approach belonging to the latter is the 
adiabatic time-dependent Hartree-Fock-Bogoliubov (AT- 
DHFB) theory started in late 1970's for the description 
of slow collective motions like low-frequency quadrupole 
vibrations and fissions, which exhibit strong non-linearity 
[Tl ] . Various versions of the ATDHFB theory have been 
proposed, e.g., by Baranger-Veneroni [l^, Villars [20l |. 
and Goeke-Reinhard However, the ATDHFB ap- 

proaches encountered some difficulties, e.g., in uniquely 
determining the collective path (see [24I for a review) , so 
that application of the theory without introducing some 
additional approximations to the real nuclear structure 
has not yet been attained. 

Still, challenge to develop a workable microscopic 
method of describing large-amplitude collective motions 
based on the ATDHF theory has been pursued. Libert 
et al.jl^ developed a practical approach assuming the 
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quadruple operators as collective coordinates and using 
the cranking mass. Quite recently, this approach was 
used in the discussion on low- lying states of 68,70,72 gg 
0. Using the generalized valley equation and the local 
RPA equation, which are based on the ATDHFB theory, 
the shape mixing in ^^Se was studied by Almehed and 
Walet [3, UHl- It is not clear, however, how the number- 
fluctuation degrees of freedom are decoupled from the 
large-amplitude shape vibrations in their work. 

On the basis of the time-dependent Hartree-Fock- 
Bogoliubov (TDHFB) theory, the self-consistent collec- 
tive coordinate (SCC) method was proposed to describe 
the large-amplitude collective motions in superconduct- 
ing nuclei [la, H^l A new scheme of solving the SCC 
equations using an expansion in terms of the collective 
momentum, called adiabatic SCC (ASCC) method, was 
formulated for describin g sh ape coexistence dynamics in 
superconducting nuclei [28l . |29|. It was firstly applied 
to the solvable multi-0(4) model to demonstrate that it 
provides an efflcient scheme to determine the collective 
path [3Qi. 

In the previous series of our works, the ASCC 
method was applied to the oblate-prolate shape coex- 
isting/mixing phenomena in ^*Se and ^^Kr, and the one- 
dimensional collective path was successfully determined 
(31I [3^ . It was shown that the triaxial deformation plays 
a crucial role in the shape mixing dynamics of these nu- 
clei. Furthermore, we constructed a four-dimensional 
collective Hamiltonian which can describe the coupled 
motion of one-dimensional collective vibration and the 
three-dimensional rotational motion of a triaxial shape. 
By requantizing the collective Hamiltonian, excitation 
spectra and quadrupole transition properties were evalu- 
ated [ssj . 

The advantage of using the ASCC method for the de- 
scription of shape coexistence dynamics is that a few col- 
lective degrees of freedom relevant to the collective mo- 
tion of interest can be extracted self-consistently from 
the TDHFB phase space. Since the collective dynamics 
is described in terms of single or a few collective vari- 
ables, it yields a clear physical interpretation of the col- 
lective dynamics. From the collective path determined 
by the ASCC method, the direction of the collective mo- 
tion can be visualized by projecting the collective path 
onto the (/3, 7) quadrupole deformation plane. It is also 
easy to evaluate the collective inertial functions (collec- 
tive mass) with respect to the 7) deformation coordi- 
nates. The obtained collective mass includes both con- 
tributions from the time-even and time-odd components 
of the moving mean-field [T8| . The time-odd contribu- 
tion from the moving mean-field is ignored in the Inglis- 
Belyaev cranking formula for the collective mass, which is 
widely used for the description of large- amplitude collec- 
tive motions. In the previous paper |34| . we have shown 
that the quadrupole pairing interaction enhances the col- 
lective mass through the time-odd component of the mov- 
ing mean-field. 

The major purpose of this paper is to give a micro- 



scopic description, on the basis of the ASCC method, of 
the oblate-prolate shape mixing dynamics in proton-rich 
Se isotopes. We show that the deformation degree of free- 
dom breaking axial symmetry plays a crucial role in the 
shape mixing. Taking into account the coupling of the 
large-amplitude shape vibrations connecting the oblate 
and prolate shapes and three-dimensional rotations of the 
triaxial shape, we show that the shape mixing gradually 
becomes weak as the rotational angular momentum in- 
creases. Dynamical reason of this trend will be clarified. 

This paper is organized as follows. In Sec. |TT1 the ba- 
sic equations of the ASCC method are summarized. In 
Sec. mil the theoretical scheme of deriving the quan- 
tum collective Hamiltonian and solving the collective 
Schrodinger equation is described. In Sec. IIVI results 
of the calculation for proton-rich 68,70,72gg isotopes are 
presented and discussed. Conclusions are given in Sec. 

El 

II. THE ASCC METHOD 

In this section the basic equations of the ASCC method 
are summarized. Details of their derivations are given in 
Ref. [H]. 

The starting point is the time-dependent variational 
principle for a TDHFB Slater determinant representing 
the collective state \(j){t)) 

5{m^§i-H\m^o, (1) 

where H denotes the microscopic Hamiltonian. In the 
SCC method, it is assumed that the collective motion 
could be described by a few canonical sets of collective 
variables. In the present application to the shape coex- 
istence phenomena, we assume that the shape dynamics 
can be described by a single collective coordinate q and 
its canonically conjugate momentum p. Since the system 
is superconducting, we also introduce the gauge angles 
(p = {(p^'^\(p^P^) together with the number fluctuation 
variables n = (n^"'',n^P') of neutrons and protons. Thus 
the TDHFB state \(l>{t)) is written in terms of these col- 
lective variables as follows. 

(2) 

where A^*^"^) = TV^"^) — A^q^'' are the number-fluctuation 

(r) 

operators about the expectation values Nq \ t denoting 
n or p. Using the generalized Thouless theorem, the in- 
trinsic state for the pairing rotation, \(j){q,p,n)) , can be 
written in terms of the moving- frame HFB state {(piq)) 
as 

|0(g,p,n))=e''^(?.f>") 10(g)), (3) 

where G{q,p,n) is a one-body operator. Note that this 
state reduces to \4>{q)) for p = and n = 0; namely. 
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\(j){q,p = 0,n = 0)) = \(f){q)). Assuming the adiabaticity 
of the large-amphtude coUective motion and the number 
fluctuations, the operator G{q,p,n) is expanded up to 
first order with respect to p and n^'^\ 



G(g,p,n)=pg((7)+^nMe(^)(g), 



(4) 



Qiq) =Q^{q) + Q''{q) 

+Qf/3(9)4«/3) , (5) 

e(^)(<z) =E (eS^(g)at at + eS^*(g)a,a„) , (6) 

a/3 

where the quasiparticle creation and annihilation opera- 
tors, and defined with respect to the moving- 
frame HFB state \<p{q)) which satisfies the vacuum con- 
ditions Qa Iff'iq)) = for them. Therefore these quasipar- 
ticle operators are also functions of the collective coordi- 
nate q. Note that the operator Q{q) contains, in addition 
to the ^-part (the first and the second terms of Eq. (O), 
the -B-part (the third term) in order to satisfy the gauge 
invariance of the ASCC equations. They are uniquely 
determined by imposing the condition [N'^'^\Q{q)] = 

The collective Hamiltonian is given by 

3 1 

n{q,p,n,I) ^{^{q,p,n)\H\^{q,p,n)) +Y,^j-^lf 
=l-(g) + ii3(q)p2+^A(-)(g)nM 



where 



1 



V{q) =n{q,p,n,I) 



p=0,n=0,f=(3' 



B{q) = 



dp^ 

dn 



p=0,n=0,/=o' 



p=0,n=0,/=o' 



(7) 

(8) 
(9) 
(10) 



are the collective potential, inverse of the collective iner- 
tial function, and the chemical potentials. The rotational 
energy term is introduced in order to treat the large- 
amplitude shape vibration and the three-dimensional ro- 
tation of triaxially deformed mean-field in a unified man- 
ner. 

The moving-frame HFB equations 

6(4>iq)\HMiq)\Hq))^0, (11) 
and the moving-frame QRPA equations 

S{c^{q)\ [HM{q),Q{q)] ~ \B{q)P{q) \<j>{q)) - 0, (12) 



S{cj>{q)\[HM{q).\p{q)]-C{q)Q{q) 

1 r r - dv ~ 

HM{q),^Q{q) 



2B{q) 



,Q{q) 



(13) 



are the basic equations which determine the collective 
path in the TDHFB phase space. They are derived by 
expanding the TDHFB equation of motion ([1]) up to sec- 
ond order with respect to p. Here Hmiq) denotes the 
moving-frame Hamiltonian 



dV 



HM{q) = H-Y^ A(^) ((7)iV(-) - —Q{q) 



The operator P{q) is defined by 

Piq)\m)^^§-^\H'i))- 

The stiffness parameter C{q) is given by 



C{q) 



1 dBdV 



dq^ 2B{q) dq dq 



(14) 



(15) 



(16) 



and connected to the moving-frame QRPA frequency as 
Lo\q) = B{q)C{q). 

The basic equations of the ASCC method are scale 
invariant, in other words, the arbitrary scale for the col- 
lective coordinate can be chosen [H]. We fix the scale 
by the condition B{q) = 1. Note also that the method 
is formulated in a gauge-invariant way; that is, the basic 
equations are invariant under the following transforma- 
tions (H. 



Q{q) ^Q(g)+aMiVW, 



(17) 



A(^H<?)^AM(<z)-aM^(g), (18) 



-{q) -^—^ — (<?) - Q-^ 'C{q). 



dq 



dq 



(19) 



Therefore, it is necessary to fix the gauge when we solve 
the ASCC equations. We adopt the same gauge fixing 
condition as in Ref . '29*1 , which is convenient to describing 
shape coexisting/mixing phenomena. 

In the following, we summarize the procedure of solv- 
ing the ASCC equations starting from one of the solu- 
tions of the static HFB equations, which corresponds 
to a local minimum of the collective potential. The 
lowest frequency QRPA eigenmode at the starting HFB 
state |(/)(g = 0)) determines the operators Q{q = 0) and 
P{q = 0). We solve the moving- frame HFB equation pT|) 
and the moving- frame QRPA equations, p2)) and (fT3|) . 
off the equilibrium to obtain the solution at q. At non- 
equilibrium HFB states, these ASCC equations are cou- 
pled with each other, so that the self-consistency between 
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the moving-frame HFB state lipiq)) and the moving- 
frame QRPA mode Q{q) is required. Let us assume that 
the sohition of the ASCC equations a,t q — 6q is already 
known. We find the solution at q by starting from solving 
the moving-frame HFB equation with the initial guess for 
the collective coordinate operator Q{q) 

Q(<Z)(") = (1 - e)Qiiq ~ 6q) + eQ2{q - Sq), (20) 

where e is a small number which mixes the lowest and 
the second-lowest solutions of the moving-frame QRPA 
equations at <; — Sq. These two solutions usually possess 
different fC-quantum numbers when the HFB mean field 
is axially symmetric. Therefore this choice for the ini- 
tial guess is crucial to find a symmetry-breaking solution 
if the previous moving- frame QRPA mode Qi{q — 5q) 
possesses the axial symmetry |32j . In this paper, we set 
£ = 0.1 in numerical calculation. 

After constructing the collective path, we evaluate the 
three rotational moments of inertia Ji (q) . For this pur- 
pose, we solve the following Thouless-Valatin equations 
at every point q on the collective path using the moving- 
frame HFB state \4i{q)) 

6{q^{q)\ [i?M(g),*.(g)] - ^Jf'm \m) = 0, (21) 

(^(g)|[vl/,(g),/,] 10(g)) (22) 

In this way, we derive the collective Hamiltonian ([7]) from 
the microscopic Hamiltonian H, which simultaneously 
describes the large-amplitude shape vibration and the 
three-dimensional rotation. 



III. REQUANTIZATION OF THE COLLECTIVE 
HAMILTONIAN 

Since the collective Hamiltonian ([7]) derived by the 
ASCC method is a classical one, it is necessary to re- 
quantize it in order to obtain collective wave functions 
describing shape mixing and discuss experimental ob- 
servables such as excitation spectra and electromagnetic 
transition probabilities. 

The total kinetic energy of the coupled motion of the 
one-dimensional large-amplitude shape vibration and the 
three-dimensional rotation is given by 

-B^^{q)q'^ + ^-Ji{q)uj^ = - ^ G'„„((7)d™d„, 

i—1 m,n 

(23) 

where LOi are angular velocities, a = (g, wi, a;2, W3), and 
the metric G„tn{q) = Smn{B-^{q), Ji{q), J2{q), Jsiq)). 
The kinetic energy term is requantized by means of the 



Pauli prescription: 

mn " 

2dq ^^'dq 4 dq \G{q)\ dq ^ 2J,{q)' 

(24) 

where \G{q)\ = B-\q)Ji{q)J2{q)J3iq) is the determi- 
nant of Gmn{q)- In this paper, we take into account the 
second term containing the derivative of \G{q)\. which 
was ignored in our previous work [s^l- Concerning the 
collective mass B~^{q), we can set it to unity without 
loss of generality, because it merely defines the scale for 
measuring the length of the collective path 28J. The 
three components li of the angular momentum operator 
are defined with respect to the principal axes (1,2,3) = 
{x',y',z') associated with the moving- frame HFB state 
(/)(q)). Care is needed when the collective path partially 
runs with axially symmetric shape where the moment of 
inertia about the symmetry axis vanishes. We discuss 
this problem in subsection IIV Cl with the concrete exam- 
ples of the collective path for ^'^Se and ^^Se. 

The collective Schrodinger equation is thus given 

(f -I- Viq))^iMk{q, = Ei^k'^iMkiq, n), (25) 

where iMk{q, ^) represents the collective wave function 
in the laboratory frame. It is a function of the collective 
coordinate q and the three Euler angles fi, and speci- 
fied by the total angular momentum /, its projection M 
on the laboratory z-axis, and the index k distinguishing 
different quantum states having the same I and M . 

Using the rotational wave functions X'^,j^(ri), the col- 
lective wave functions in the laboratory frame is written 
as 

I 

^iMk{q, ^)^Y. ^iKk{q)mMK), (26) 

K=0 



"V l67r^(l + ^Ko) ^''^'^^^"^ + (-)X^-k(^^)), 

(27) 

where ^iKk{q) represents the large-amplitude vibra- 
tional motion, and the sum in Eq. (j26p is restricted to 
even K. This form of the collective wave function is 
analogous to that of the Wilets-Jean 7-unstablc model 
but the role of the triaxial deformation variable 7 
is played here by the collective coordinate q. Needless to 
say, this is a particular form in the general framework of 
the Bohr and Mottelson [3^. 

Normalization of the vibrational part of the collective 
wave functions is given by 

/ dr' ^*iKk {q)^iKk' {q) = Skk' , (28) 
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where the volume element is 

dr = dr'dn = ^/\G{q)\dqdn. (29) 

The boundary conditions for the collective Schrodinger 
equation (pS)) can be specified by projecting the collective 
path obtained by the ASCC method onto the 7) plane 
and by using the well-known symmetry properties of the 
Bohr-Mottelson collective Hamiltonian j33|, i3^, i37|] . 

IV. SHAPE MIXING IN PROTON-RICH SE 
ISOTOPES 

A. Model Hamiltonian and parameters 

For the microscopic Hamiltonian, we use a version of 
the pairing-plus-quadrupole (P+Q) force model [H, [3§| 
which includes the quadrupole-pairing in addition to the 
monopole-pairing interaction. Two major shells (A^sh = 
3, 4) are considered as the active model space for neu- 
trons and protons. The single-particle energies are cal- 
culated using the modified oscillator potential [i^. As 

in Ref. [s^, the monopole-pairing strength G^J^ and the 
quadrupole-particle-hole interaction strength x' for ^^Se 
are determined such that the magnitudes of quadrupole 
deformations and pairing gaps at the oblate and prolate 
local minima approximately reproduce those obtained in 
the Skyrme-HFB calculation by Yamagami et al. 
The interaction strengths for ™^^^Se are then determined 
from those of ^*Se, assuming a simple mass number de- 
pendence G^"") - and x' For the 
quadrupole-pairing interaction strengths g'^\ the self- 
consistent values proposed by Sakamoto and Kishimoto 
[4l| are evaluated from the monopole pairing interaction 
(see [ll] for details). These values of the interaction 
strengths are listed in Table HI 

Following the conventional treatment of the P+Q 
model, we ignore the Fock term, so that, in the follow- 
ing, we use the abbreviation HB (Hartree-Bogoliubov) in 
place of HFB. The effective charges (e„, ep) — (0.4e, 1.4e) 
are used in the calculation of E2 transition matrix ele- 
ments. In numerical calculation of solving the ASCC 
equations, we use Sq = 0.01. 

B. Properties of the HB states and the QRPA 
vibrations 

The static HB solution and the QRPA calculation 
based on it provide the ASCC solution a.t q — Q. Prop- 
erties of the HB mean-field and of the QRPA modes are 
summarized in Table [Til In all the three isotopes, we ob- 
tain two HB solutions possessing the oblate and prolate 
shapes. While the magnitude of the quadrupole defor- 
mation of the oblate HB state depends on the neutron 
number rather weakly, that of the prolate HB state signif- 
icantly increases from ^*Se to ^^Se. This trend of equi- 
librium deformation is consistent with what we expect 



TABLE L Strengths of the monopole-pairing, the quadrupole 
particle-hole, and the quadrupole-pairing interactions 
adopted in the numerical calculation. The same monopole- 
pairing strength is used for neutrons and protons, i.e., we set 
Go = Gq"'' = Gq'\ For the quadrupole-pairing interactions, 
the strengths multiplied by the oscillator length biquadrate, 

Ga^^' = Gi^^h'^, are shown. 



Go (MeV) x' (MeV) G^'"' (MeV) g'J-"^ (MeV) 



-^^Se 


0.320 


0.248 


0.185 


0.185 


™Se 


0.311 


0.236 


0.174 


0.184 


^^Se 


0.302 


0.225 


0.161 


0.183 



from the deformed shell gap in the Nilsson diagram. The 
oblate HB solutions are always the lowest in energy, but 
the energy difference between the oblate and prolate HB 
local minima are only 0.3 ~ 0.6 MeV. 

Concerning the QRPA vibrations in ^^Se, the 7- 
vibration is the lowest frequency mode, and the /3- 
vibration is the second-lowest mode both at the oblate 
and prolate minima. The situation is different in ™Se, 
where the /^-vibration is the lowest mode both at the 
oblate and prolate minima. In the case of ''^Se, the 
7-vibration is the lowest mode at the oblate minimum, 
while the /3-vibration is the lowest mode at the prolate 
minimum. 

It is also seen in Table Hi] that the rotational moments 
of inertia at the prolate minimum significantly increases 
from ^*Se to ^^Se following the increase of the quadrupole 
deformation /3. The calculated values for ^^Se and ™Se 
seem a little too small, however, compared to the values 
experimental data suggest. We plan to make a more 
detailed analysis about this problem in future. 



C. Collective path 

We have solved the ASCC equations and determined 
the collective path choosing one of the HB solutions in 
Table HIl in each nucleus and setting it as \(f){q — 0)). The 
results are displayed in Fig. [1] where the obtained collec- 
tive paths are drawn in the {(3, 7) deformation plane. 



In this nucleus, the potential barrier height is about 0.5 
and 0.07 MeV measured from the oblate and prolate local 
minima, respectively. The oblate HB state is chosen as a 
starting state for solving the ASCC equations. Though 
the collective path for ^^Se is already reported in the 
previous work [33| . we summarize the character of the 
solution of the ASCC equations for later convenience. 
The collective path starting from the oblate HB states 
almost follows the triaxial potential valley. 
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TABLE II: Calculated values for the quadrupole deformations (/3,7), the monopole pairing gaps for neutrons and protons 
(A("\ Aq'''), the potential energy V measured from the lowest minimum of the HB equilibrium states, the frequencies {lUj,uji3) 
of the lowest two QRPA modes at the HB equilibrium states, the collective mass M for the lowest QRPA mode. The QRPA 
modes with K = and K — 2 are denoted (3 and 7, respectively, where K is the symmetry axis component of the vibrational 
angular momentum. The rotational moments of inertia J about the axis perpendicular to the symmetry axis are also shown. 



/3 7 AIP^ (MeV) A^,''^ (MeV) V (MeV) uj^ (MeV) top (MeV) M (MeV-^) J (MeV"') 



""Se (ob) 0.296 60° 


1.17 


1.26 





1.373 


2.131 


50.96 


6.38 


^**Se (pro) 0.260 0° 


1.34 


1.40 


0.41 


0.886 


1.367 


34.29 


4.60 


■'"Se (ob) 0.313 60° 


1.21 


1.16 





1.617 


1.421 


83.07 


7.52 


™Se (pro) 0.325 0° 


1.34 


1.30 


0.55 


1.161 


1.120 


47.51 


6.89 


''^Se (ob) 0.268 60° 


1.42 


1.16 





1.294 


1.482 


52.90 


6.18 


■^^Se (pro) 0.381 0° 


1.08 


1.23 


0.32 


1.411 


1.042 


72.28 


10.25 





pcosy (Jcosy Pcosy 

FIG. 1: (Color online) The collective paths for ®*~^^Se obtaind by the ASCC method. The collective path projected onto 
the (/3, 7) deformation plane are drawn by solid lines on the potential energy surface. The equipotential lines are drawn every 
100 keV. Note that the collective path is symmetric with respect to the reflections about the prolate (7 = 0°) and the oblate 
(7 = 60°) axes. The collective path going along the symmetry axis eventually terminates at large /3 when the neutron pairing 
gap collapses (see the text). 



In Fig. [2I we see that the deformation /3 almost stays 
constant during when the triaxial deformation 7 changes 
from 60° to 0° along the collective path. This clearly 
indicates that the triaxial degree of freedom plays much 
more important role than the axial degree of freedom in 
^*Se. The 7-dependence of the calculated moments of 
inertia exhibits a behaviour similar to the irrotational 
moments of inertia; two of them coincide at the axially 
symmetric limit while the largest moment of inertia is 
about the intermediate axis. 

In Fig. [21 the collective mass defined as a function of 
the geometrical length ds = \J dp^ + P^d'y^ in the (/3, 7) 
plane, 

M{s{q)) = B-\q) I (f )' + (^)'} , (30) 

is also presented. As mentioned in Sec. IIIIl we can set 
M{q) = B{q)-^ = 1 MeV"^ using the units where h=l 
and the collective variables {q-,p) are non-dimensional. 



In this nucleus, the potential barrier height is about 
0.7 and 0.1 MeV measured from the oblate and prolate 



local minima, respectively. The collective path is ob- 
tained starting from the prolate HB state. The two HB 
local minima are connected by the one-dimensional path. 
Since the QRPA mode with the lowest frequency at the 
prolate shape has /3- vibrational character with K = 0, 
the collective path starting from the prolate HB state 
goes along the axial symmetry axis in the beginning. As 
seen in Figs. [T] and [3l at q ~ 0.4 (/3 ~ 0.27), the collective 
path deviates from the 7 = 0° line due to the character 
change of the lowest mode from K = to K — 2. To 
describe such a dynamical breaking of the axial symme- 
try taking place along the collective path, it is crucial to 
use Eq. (PO)) as an initial trial for self-consistently deter- 
mining the collective coordinate operator Q{q). The col- 
lective path encounters a similar avoided crossing of the 
moving- frame QRPA modes at q ~ 1.8 (the oblate side 
with f3 ~ 0.27). Then, the triaxial path again changes its 
direction and go along the 7 = 60° line. This kind of dy- 
namical breaking of the axial symmetry was previously 
reported in the analysis of the collective path for ''^Kr 
[33 . [ssj . At the oblate side of the collective path, the /3- 
vibrational degrees of freedom strongly couples with the 
pairing-vibrational degree of freedom of neutrons, and it 
ends at a large /3 point where the neutron pairing gap 
collapses. When approaching this point, the collective 
mass diverges. 

As the rotation about the symmetry axis disappear. 
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FIG. 2: Results of the ASCC calculation for ®*Se. The monopole pairing gaps Aq (g), the quadrupole pairing gaps A^o {q) 
and A22''(fj'), the collective potential V{q), the collective mass M{s{q)), the rotational moments of inertia Ji{q), the lowest two 
moving- frame QRPA frequencies squared a;^(g), the axial quadrupole deformation and the canonical collective coordinate 

q are plotted as functions of 7(5). 



the moments of inertia Jz{(l) and ^^2(9) should be 
dropped in the determinant |G(g)| of the metric Gmnil) 
when the collective path runs along the 7 = and 
7 = 60 lines, respectively. To avoid discontinuity at 
the point where the collective path starts to deviate from 
the symmetry axis, we use \G{q)\ = B~^{q)Ji{q) for the 
whole region of the collective path of this nucleus. 

In the region of prolate shape with f3 > 0.34, the lowest 
two QRPA modes with f3 and 7 characters approximately 
degenerate in energy and compete (see Fig. [3]) . In such a 
situation, it may be appropriate to introduce two collec- 
tive coordinates. In the present calculation, however, we 
have solved the ASCC equations in this region assuming 
that the collective path continues to go along the 7 = 0° 
axis. We shall attempt to introduce two collective coor- 
dinates in the same framework of the ASCC method in 
future. 



and 01 the collective path directs to the triaxial region 
because the character of the lowest QRPA mode at the 
oblate minimum is 7-vibrational. At g ~ 0.2 in the triax- 
ial region, the collective path curves due to an interplay of 
the lowest two moving-frame QRPA modes. The collec- 
tive path reaches the prolate side at q ~ 1.6 (/3 ~ 0.32). 
Then the path changes its direction to the 7 = 0° line. 
Thus, the oblate and prolate local minima are connected 
by a single collective coordinate. After passing through 
the prolate minimum at g ~ 2.1 (/3 ~ 0.38), it contin- 
ues to go along the 7 = 0° line and finally terminates at 
f3 ~ 0.42 where the neutron pairing gap collapses. Corre- 
spondingly, the collective mass increases with increasing 
/3 and finally diverges. 



In this nucleus, the potential barrier height is about 0.5 
and 0.3 MeV measured from the oblate and prolate local 
minima, respectively. The collective path is determined 
starting from the oblate HB state. As seen in Figs. [1] 



As discussed above for '''^Se, the moment of inertia 
^73(9) should be dropped in the determinant |G(9)| when 
the collective path runs along the 7 = line. Since the 
collective path for ^^Se does not run along the 7 = 60 
line at all, we use \G{q)\ = B~^{q)Ji{q)J2{q) for the 
whole region of the collective path of this nucleus. 
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FIG. 3: The same as Fig. [2]but for '^'^Se plotted as functions of q along the collective path. The point g = corresponds to the 
prolate local minimum, while the oblate minimum is located at g = 2.18. These positions are indicated by arrows with P or O. 



D. Shape mixing, excitation spectra, quadrupole 
transitions and moments 



We have calculated collective wave functions solving 
the collective Schrodinger equation (^5)1 and evaluated 
excitation spectra, quadrupole transition probabilities 
and spectroscopic quadrupole moments. Below we dis- 
cuss these results denoting the eigenstates belonging 
to the ground and excited bands as 0^,2^,4^,6^ and 
, 2 J , 4+ , 6 J , respectively. 



1. 



In Fig. [51 excitation spectrum and B(E2) values calcu- 
lated for ^^Se are displayed together with experimental 
data. It is seen that the calculation yields two bands 
which exhibit significant deviations from the regular ro- 
tational pattern. We note, in particular, that the cal- 
culated state is located above the 2^ state. This 
is consistent with the available experimental data where 
the state has not yet been found. Wc see in Fig. [HI 
that the vibrational wave functions of the Ol", Oj, 2f and 
2^ states spread over the whole extent of 7 from the 
oblate to the prolate shapes. This result of calculation 
is reasonable considering the very low potential barrier 
along the triaxial collective path, as we have already seen 



in Fig. [21 The unusual behavior of the excited 0+ state 
noted above suggests that the low- lying states in ^*Se are 
in an intermediate situation between the oblate-prolate 
shape coexistence and the Wilets-Jean 7-unstable model 
[ssj . In fact, we can find a pattern in the calculated E2- 
transition probabilities, which is characteristic to the 7- 
unstable situation; for instance, -B(E2; 6^ — > 6^), B(E2; 
4 J 4f ) and B{E2; 2 J ^ 2f) are much larger than 
B{E2; 6j ^ 4+), B{E2; 4+ ^ 2+) and B{E2; 2+ ^ 0+). 
This point will be discussed with a more general perspec- 
tive in a future publication [43| . It is quite interesting to 
notice that the shape mixing becomes weak as the an- 
gular momentum increases, and the collective wave func- 
tions of the 4+ and 6+ states tend to localize in the region 
near either the oblate or the prolate shape. Namely, it 
becomes possible to characterize the 4+ and 6'^ states as 
oblate-like or prolate-like. 



>Se 



Calculated and experimental excitation spectra and 
B(E2) values for ™Se are displayed in Fig. [71 The ex- 
citation energies of the ground band is well reproduced. 
The calculated B(E2;2+ ^ Of) value 390 fm^ is also 
in reasonable agreement with the experimental data 342 
fm''. The calculated E2-transition probabilities ex- 
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FIG. 4: The same as Fig.[2]but for ^'^Se plotted as functions of q along the collective path. The point g = corresponds to the 
oblate minimum, while the prolate local minimum is located at g = 2.34. These positions are indicated by arrows with O or P. 
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FIG. 5: Theoretical (/e/i) and experimental (right) excitation spectra and B{E2) values for low-lyine states in Se. Only 
_B(E2)'s larger than 1 Weisskopf unit are shown in units of fm''. Experimental data are taken from [y, In- 



hibit a pattern somewhat different from that of ^^Se; 
for instance, we see significant cross talks among the 
2]^, 2^, 4^ and 4 J states. The vibrational wave func- 
tions of the of , Oj, 2f and 2 J states displayed in Fig. [8] 
show strong oblate-prolate shape mixings. In contrast, 
the 4 and 6^ (42 and 6 J) states are rather well local- 



ized around the prolate (oblate) shape. Thus, the char- 
acteristic cross talks of the E2 transition strengths men- 
tioned above are associated with the significant change in 
localization properties of the vibrational wave functions 
between angular momenta 2 and 4. Experimental data 
for such inter- and intra-band i?(E2) values will certainly 
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oblate minimum, the moment of inertia J ~ 10.3 MeV~^ 
at the former is appreciably larger than J7 ~ 6.2 MeV"'^ 
at the latter. The difference between the rotational ener- 
gies at the prolate and oblate minima is easily evaluated 
to be about 0.19, 0.64, 1.35 MeV for the 2+, 4+, 6+ states, 
respectively. Therefore, the prolate shape is favoured in 
order to reduce the rotational energy. As the rotational 
angular momentum increases, this rotational effect be- 
comes more important and overcomes the small disad- 
vantage in the potential energy. In contrast, for the 0+ 
ground state where the rotational effect is absent, the vi- 
brational wave function takes the maximum at the oblate 
minimum. For the 2\ state, the difference of the rota- 
tional energies, about 0.19 MeV, is slightly smaller than 
that of the potential energies, about 0.32 MeV, so that 
its collective wave function exhibits a transitional char- 
acter from oblate-like to prolate-like. For the excited 
states, the vibrational wave functions possess the domi- 
nant bumps around the oblate shape, exhibiting at the 
same time the second bumps around the prolate shape. 



FIG. 6: (Color online) Vibrational wave functions ^iKk{q) 
squared of the lowest (left) and the second-lowest states (right) 
for each angular momentum in ®*Se. In each panel, different 
K components of the vibrational wave function and the sum 
of them are plotted as functions of 7(g). 



]uadrupole moments 



serve as a very good indicator of the shape mixing. 
3. '^Se 

Calculated excitation spectrum and _B(E2) values for 
^^Se are shown in Fig. [Altogether with experimental data. 
It is seen that the experimental spectrum is reproduced 
fairly weU. The calculated B(E2;2+ ^ 0+) value 390 
fm"* is also in good agreement with the experimental data 
405 fm**. We see that the calculated inter-band E2 
transitions, B(E2; 4+ ^ 4+), B(E2; 4+ ^ 2+), B(E2; 
2\ 2^) and i3(E2; 4+ ^ 2^), are reduced from those 
in ™Se, except for B(E2; 0+ 2+). 

The vibrational wave functions are displayed in Fig. 1101 
Similarly to "^"Se, the 0'^ wave function widely spreads 
over the triaxial region. It takes the maximum at the 
oblate shape, but extends to the prolate region. The 2^ 
wave function also extends the whole region of 7. In 
the ground band, the prolate character develops with in- 
creasing angular momentum, as clearly seen in the wave 
functions of the and 6^ states. 

Dynamical reason why the prolate character of the 
ground band develop with increasing angular momentum 
may be understood in terms of the competition between 
the potential and kinetic energies as function of the col- 
lective coordinate q. We find that the rotational energy 
term plays a particularly important role. Because the 
quadrupole deformation /3 ~ 0.38 at the prolate local 
minimum is much larger than that (/3 ~ 0.27) at the 



The spectroscopic quadrupole moments calculated for 
68,70,72sg j^j.g displayed in Fig. [IH For ^^Se, the 4+ and 
6^ states possess positive signs indicating dominance of 
the oblate character, while the 4^ and 6^ states have 
negative signs indicating dominance of the prolate char- 
acter. In contrast to the ^^Se case, the 4^ and 6^ states 
in ''"Se and ^^Se have negative signs indicating the growth 
of prolate character of these states with increasing rota- 
tional angular momentum. These results are in qualita- 
tive agreement with the HFB-based configuration-mixing 
calculation reported by Ljungvall et al. Both calcu- 
lations indicate the oblate (prolate) dominance for the 
ground (excited) band in ^'^Se while the prolate charac- 
ter develops with increasing angular momentum for the 
ground bands in ™Se and ''^Se. For the excited bands 
of ^°Se and ^^Se, results of the configuration-mixing cal- 
culation are not reported in Ref. [g]. Our calculation 
indicates the growth of oblate character for the the 4^ 
and 62 states in these isotopes. 

Careful interpretation is necessary when absolute val- 
ues of the calculated spectroscopic quadrupole moment 
are small. In our results of calculation, small values have 
nothing to do with spherical character of the states of 
interest; it is a particular consequence of large-amplitude 
shape vibration. Namely, we find a number of situations 
where the contributions from the components of the vi- 
brational wave function with 7 > 30° are largely can- 
celed with those from 7 < 30° . Such a cancellation is the 
main reason why the calculated quadrupole moments are 
rather small for all the 2^ and 2j states of interest. 
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FIG. 8: (Color online) Same as Fig. |6] but for ™Se. The 
vibrational collective wave functions squared are plotted as 
functions of q. The arrows indicate the positions of the oblate 
(O) and the prolate (P) minima. 



states for ^^Se from that in Ref. [s^. This indicates im- 
portance of proper treatment of the volume element. In 
numerical calculations for ^°Se and "^^Se, however, the 
volume element was treated in an approximate way. We 
plan to examine the validity of this approximation by 
deriving a five-dimensional quadrupole collective Hamil- 
tonian on the basis of the ASCC method and make a 
detailed comparison of the present results with those of 
the five-dimensional calculation [4^ . 

Another question is the validity of evaluating the ro- 
tational moment of inertia after determining the collec- 
tive path, Obviously, the assumption that the collective 
path does not change due to the rotational motion will 
be eventually violated with increasing angular momen- 
tum. Namely, the present approach may be valid only for 
low-spin states. We have therefore restricted our calcu- 
lation to low-spin states with angular momentum / < 6. 
By comparing with the five-dimensional calculation men- 
tioned above, we shall be able to examine also the range 
of applicability of the present approach. Alternatively, 
one can use the rotating mean field when determining 
the collective path. Such an approach was once tried in 
Refs. m and in. 



V. CONCLUSIONS 



5. Discussion 

Before concluding, we remark on a few questions to be 
examined in a future publication. 

In this paper, we have taken into account the function 
G{q) in the volume element while it was put unity 
in the previous calculation [33j. Thus, we have obtained, 
for instance, different ordering between the 0^ and 2j 



Using the ASCC method we have investigated the 
oblate-prolate shape coexisting/mixing phenomena in 
proton-rich selenium isotopes, ^^'^'''^^Se. The collective 
paths connecting the oblate and prolate HB local min- 
imum were successfully determined. Requantizing the 
collective Hamiltonian obtained by means of the ASCC 
method, we have derived the quantum collective Hamil- 
tonian which describes the large-amplitude shape vibra- 
tion along the collective path and the three-dimensional 




FIG. 9: The same as Fig. [5] but for ^^Se. Experimental data are taken from [^.14^. 
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FIG. 10: (Color online) Same as Fig. H but for ^^Se. The 
vibrational wave functions squared are plotted as functions of 
q. The arrows indicate the positions of the oblate (O) and 
the prolate (P) minima. 



rotational motion in a unified manner. Solving the col- 
lective Schrodinger equation, we have calculated exci- 
tation spectra, E2 transition probabilities and spectro- 
scopic quadrupole moments. It has been shown that the 
basic properties of the coexisting two rotational bands in 
low-lying states of these nuclei are well reproduced. 

The result of calculation clearly shows that the oblate- 
prolate shape mixing becomes weak as the rotational an- 
gular momentum increases. We have analyzed dynamical 
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FIG. 11; Spectroscopic quadrupole moments calculated for 
the low-lying states in ^^'^"'^^Se. The square, circle and trian- 
gular symbols represents those for ^*Se, ^"Se and ^^Se, respec- 
tively. The filled (open) symbols denotes the ground (excited) 
band. 



origin of this trend and found that the rotational energy 
plays a crucial role in determining the degree of localiza- 
tion of the collective wave function in the (/3, 7) deforma- 
tion space. The rotational effect causing the localization 
of the collective wave function may be called "rotational 
hindrance of shape mixing." To our knowledge, impor- 
tance of such a dynamical effect has not been received 
enough attention in connection with the oblate-prolate 
shape coexistence phenomena of interest. The rotational 
hindrance effect will be discussed with a more general 
perspective in a future publication (4^ . 
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